Exchange in silicon-based quantum computer architectures 
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The silicon-based quantum computer proposal has been one of the actively pursued ideas during 
the past three years. Here we calculate the donor electron exchange in silicon and germanium, and 
demonstrate an atomic-scale challenge for quantum computing in Si (and Ge), as the six (four) 
conduction band minima in Si (Ge) lead to inter-valley electronic interference, generating strong 
oscillations in the exchange splitting of two-donor two-electron states. Donor positioning with atomic 
scale precision within the unit cell thus becomes a decisive factor in determining the strength of the 
exchange coupling — a fundamental ingredient for two-qubit operations in a silicon-based quantum 
computer. 
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Following the seminal proposal [y by Kane there has 
been a great deal of activity in efforts to develop a 
silicon-based quantum computer (QC) architecture. The 
basic ideas of the Kane proposal are simple and attrac- 
tive: to use donor nuclear spins as quantum bits (qubits), 
and to utilize the vast infrastructure and technology as- 
sociated with the Si industry to fabricate precisely con- 
trolled Si nanostructures, where exchange effects between 
electrons and nuclei in neighboring donor impurities (e.g. 
31 P in Si) could serve as the two-qubit gates, similar to 
the electron-spin-based QC proposal by Loss and DiVin- 
cenzo ||. The motivation for a Si quantum computer is 
obvious: Once the basic one-qubit and two-qubit opera- 
tions have been demonstrated using donor impurities in 
Si nanostructures, computer chip fabrication technology 
associated with the existing and dominant Si industry 
will easily enable the scale-up of information processing 
involving large number of donor nuclear spin qubits. In- 
deed, one of the formidable stumbling blocks in develop- 
ing working quantum computer hardwares has been the 
scale-up problem, as the demonstrated qubits in trapped 
ion and liquid state NMR techniques are not readily scal- 
able in any significant manner 

A great deal of experimental work is currently be- 
ing aimed at developing suitable qubits in Si nanostruc- 
tures with precisely introduced dopant impurities, using 
both a "top-down" approach with ion-implantation, and 
a "bottom-up" approach with MBE growth and scanning 
tunneling microscopy f§. In the Si QC model donor 
electrons act as shuttles between different nuclear spins. 
For two-qubit operations, which are required for a uni- 
versal QC, both electron-electron exchange and electron- 
nucleus hyperfine interaction need to be precisely con- 
trolled. These are unquestionably formidable experimen- 
tal problems. In the original proposal, Kane used the 
Herring-Flicker exchange formula for two hydrogenic 
centers to obtain an order of magnitude estimate of the 
electron exchange among donors in Si Q|. However, as he 
also pointed out, donor exchange in Si is not hydrogenic. 



In this Letter we show that exchange effects in pro- 
posed donor nuclei based Si QC architectures are actually 
very subtle due to quantum interference effects inherent 
in the complicated Si band structure. In particular, spe- 
cial care, going far beyond what is being currently at- 
tempted in the fabrication of Si QC, will be required in 
controlling the surface gates crucial to QC operations. 
We also provide results for Ge which seems to have cer- 
tain advantages over Si as a candidate for dopant ex- 
change based quantum computation. Unfortunately, Ge 
is not such an attractive practical alternative since there 
is no well-established Ge technology to take advantage 
of. 

The study of shallow impurities in Si and Ge is a 
quite mature field However, complications aris- 

ing from the anisotropy of the electron effective mass 
and from interference among the Bloch wavefunctions at 
the degenerate conduction-band edges were never fully 
explored in detail. Both effects were discussed by An- 
dres et al statistically, in a study of magnetic sus- 
ceptibility of Si:P system performed in the context of 
localization and magnetic phase transitions. Silicon con- 
duction band has six minima located along the <100> 
directions, at about 85% between the center (r point) 
and the boundary (X points) of the Brillouin zone (BZ): 
k± z = 0.85(0, 0, ±27r/a), and the equivalent x,y direc- 
tions. The interference between these valleys causes fast 
oscillations (on the scale of the inter-atomic spacing, 
quite different from the slowly varying exchange interac- 
tion in the spin-based quantum dot QC model P,pl]-^3[) 
in the exchange interaction. Such oscillations cannot 
be accounted for by a simple calculation using hydro- 
genic centers. Here we perform a Hcitlcr-London calcu- 
lation of the exchange coupling Jll[] between two sub- 
stitutional donors in bulk Si, incorporating the effects 
of both the anisotropic effective mass and the valley in- 
terference. Within the envelope function approach, the 
ground single valence donor electron state (with A\ sym- 
metry) is written as a symmetric superposition of the six 
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conduction band minima 
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where 4>u (r) = u A1 (r)e lk » i ' r are Bloch wavefunctions. Here 
refer to the six conduction band minima, and F^iv) 
are the corresponding envelope functions. We use the 
Kohn-Luttinger variational form for these anisotropic en- 
velopes HH, e.g., 
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To determine the effective Bohr radii a and b, we have 
used the most recently measured Si and Ge parame- 
ters |l4| ] (especially longitudinal and transverse effective 
masses) to perform a variational calculation || . The cal- 
culated wavefunction widths, as shown in Table |, are 
quite close to those obtained over 40 years ago ||. The 
large difference in transverse mass has only minimal ef- 
fect on our results. 

To calculate the exchange splitting between the ground 
singlet and triplet states for an impurity pair in Si, we 
use an approximate form of the Heitler-London approach 
[pToh , which leads to 
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where R is the relative position of the impurity nuclei 
pair, and 
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The second summation in (||) comes from the reciprocal 
lattice expansion of the periodic part of the Bloch func- 
tion, u M (r) = X^K c K ejK r ' an d is identically 1 when R 
is an fee crystal lattice vector. Integrals in Eq. (||) with 
rapidly oscillating integrands are neglected in Eq. 
For n = v, we calculate the integrals in (|^) by replacing 
the denominator |ri — r2 1 by its value for the line along 

TABLE I. Experimental values for the lattice parameter, 
dielectric constant and electron effective masses (in units of 
the free electron mass) for Si and Ge [[HJ. The last two 
columns give the calculated values of the effective Bohr radii. 
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the impurities (To) , and assume = yj j^juu when /i ^ 
v. All are excellent approximations for large separations 
(|R| > a, 6> a). 

Germanium conduction band has four minima along 
the <111> directions, at the zone boundary L points. 
Straightforward changes in (|l|) lead to expressions for 
the exchange coupling in Ge similar to (Q), with fi and 
v along the <111> directions. Because the conduction 
band minima are located at the zone boundary, the oscil- 
lation in the exchange coupling should display a simpler 
behavior than in Si when both donors are on the same 
fee sublattice, thus making Ge an intriguing candidate 
for the purpose of quantum computing. 

We calculate the exchange energy for a pair of donors 
at R^ and as a function of their relative position 
R = Ra — Rb from Eq.Q, taking the expression in the 
square bracket equal to unity. This means that the re- 
sults are appropriate for any lattice vector R, while for 
general inter-donor distances our values should be taken 
as an estimate. For deffniteness, we consider R along 
high-symmetry directions in the crystal lattice. Fig. [j] 
shows the calculated values of J(R) for donors in Si and 
Ge with R along the [100], [110] and [111] directions 
(frames (a), (b) and (c) respectively). The solid lines 
give the results for Si, and exhibit the expected decay of 
J with increasing |R| due to the decrease of the donor 
wavefunction overlap. Other general features of Si and 
Ge exchange are clearly illustrated in this figure, namely 
the oscillatory behavior of J superimposed on its overall 
decay with distance, and the strong anisotropy of J(R) 
which is apparent by comparing different frames. Both 
features are consequences of the host material band struc- 
ture and have not been considered in detail in previous 
studies pi, either for simplicity |l|,||, or because such 
effects are averaged out for a random donor distribution 
pOfl . The filled circles in Fig. [l| indicate the accessible 
values of J(R) when the impurity pair in Si is located 
at lattice sites along the considered directions. Since for 
Si the conduction band minima correspond to points in- 
side the BZ, the period of oscillation in J and the lattice 
periodicity are not commensurate. We have also inves- 
tigated the effect of small perturbations in the atomic 
positions: The open squares give the resulting exchange 
values when one of the impurities is slightly displaced 
to off-lattice positions. The set of squares around each 
circle in Fig. ^ corresponds to displacements along differ- 
ent directions, with the distance from the original lattice 
position arbitrarily taken as S = 0.235 A (i.e., 10% of 
the nearest-neighbor distance in Si). The squares fol- 
low to a very good approximation the behavior of the 
calculated J(R), with R along the unperturbed crystal 
direction, regardless of the direction of 5. The upper 
bound of AJ/J calculated from the small-displacement 
data is about 2.5 to 5% along the [100] direction, which 
means that the requisite control over donor positioning 
should be much better than 10% of the nearest neighbor 
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distance, i.e. better than 0.235 A, a rather difficult task. 

The corresponding results for donor pairs in Ge are 
given by the dashed lines, diamonds and crosses in Fig. [|. 
In these calculations we have also assumed nearly free 
electrons, which means Cg 3> c^^ . The main qualita- 
tive difference between the calculated exchange coupling 
for donors in Ge and Si arises from the different locations 
of the conduction band minima in the BZ. Since for Ge 
the minima occur at the zone boundary L points, the 
oscillations in J (JR.) are commensurate with the lattice 
periodicity. For donors at lattice sites along the [100] di- 
rection, the accessible values of J(R) correspond exactly 
to successive local maxima of Eq.(^). As a consequence, 
off-lattice displacements from the original lattice sites by 
10% of the nearest-neighbors distance in Ge (S = 0.245 
A) have negligible effect on the exchange coupling for 
donors along this direction. Quantitatively, the longer 
range for the interactions in Ge as compared to Si is due 
to the larger values of the effective Bohr radii for Ge (see 
Table |). 
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FIG. 1. Calculated exchange coupling between two phos- 
phorus donors in Si (solid lines) and Ge (dashed lines) along 
high-symmetry directions for the diamond structure (see in- 
set in Pig.^j) . Values appropriate for impurities at substitu- 
tional sites are given by the circles (Si) and diamonds (Ge). 
Off-lattice displacements by 10% of the nearest-neighbor dis- 
tance lead to the perturbed values indicated by the squares 
(Si) and crosses (Ge). 

Our results indicate that a 31 P donor array along the 
[100] direction in either of the host materials meets the re- 
quirements for quantum computer implementation. The 
exchange energy for two 31 P substitutional impurities 100 
to 200 Aapart along this direction ranges between a few 



tenth of one meV to one /ueV, corresponding to less than 
100 picosecond to 100 nanosecond gate operation time. 
These values are robust with respect to small off-lattice 
displacements. 
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FIG. 2. Variations in the exchange coupling between two 
phosphorus donors in Si (solid symbols) and Ge (open sym- 
bols). Lattice displacements into one of the first (red), second 
(green), or third (blue) neighboring positions follow the color 
code represented in the inset, which depicts a diamond struc- 
ture with colored bonds. The numbers in the parenthesis next 
to the data points are their degeneracies, respectively. For 
clarity, not all bonds connecting neighbors are shown in the 
inset. Black circles refer to the reference positions Ro, with 
the two donors exactly along the [100] direction. The red 
squares give the exchange constant when one of the donors 
is displaced into one of its nearest neighbor sites, while green 
diamonds represent the second nearest neighbor values and 
blue circles the third nearest neighbor values. The points 
connected by a line refer to pairs along the [100] direction, 
displaced by ±a with respect to the reference position, and 
may be easily identified in Fig. 0(a). 

Displacements of the 31 P atoms into neighboring lat- 
tice sites, which are bound to occur in the fabrication 
process due to either uncontrolled implantation or sur- 
face diffusion during MBE growth, also deserve careful 
theoretical investigation. In Fig. || we present the calcu- 
lated J(R) for a specific relative position Ro along the 
[100] direction (black symbols), as well as for values of 
R = Ro + A, with A ranging over the 4 nearest neigh- 
bors, 12 second nearest neighbors and 6 third nearest 
neighbors in the diamond structure (see inset of Fig. |^). 
These displacements cause a relatively small change in 
the inter-donor distance (see horizontal scale in the fig- 
ure). For third neighbors (blue symbols) the trends in 
Fig. ^ are reasonably preserved. Surprisingly, most first 
and second neighbors displacements essentially destroy 
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the exchange coupling. This effect is entirely due to the 
cos^^ • R) factors in Eq. (|J), coming from the plane- 
wave part of the Bloch functions, and is not directly re- 
lated to the numerical values or approximations involved 
in obtaining the coefficients. Therefore we do not at- 
tempt to refine the envelope function calculation at this 
stage, e.g. by including central cell corrections M and 
inter- valley coupling |n|. The results obtained here are 
inconsistent with the commonly accepted idea that when 
the Bohr radii of the envelope functions are much bigger 
than the lattice parameter, the substitutional position- 
ing of the impurity atoms is not of much importance in 
determining the exchange coupling. 

In a Si QC architecture there are surface gates (so 
called A and J gates) that are used to control single 
and two-qubit operations 0. By adjusting the gate po- 
tentials), one can shift the center of a donor electron 
wavefunction, thus quickly traversing the fast oscillating 
terrain in the exchange coupling. The results obtained 
in this paper also indicate that donor electron exchange 
depends very sensitively on the applied gate potential(s). 
Furthermore, whether the electrons are trapped by donor 
ionic potentials or by the applied gate potentials, the fast 
oscillation in exchange will persist as long as the electrons 
are in the bulk of a Si crystal. For a 2-qubit gate of a 
QC, it is the time integral of the exchange constant J 
that determines the gating time ||. The implication of 
the oscillatory exchange is that the A and J gate volt- 
ages corresponding to the peak exchange coupling have 
to be well-controlled, optimally close to a local maximum 
where the exchange is least sensitive to the gate voltage. 
Since the oscillatory exchange period is close to lattice 
spacing, the positioning of the donor electrons by the A 
and J surface gates must be controlled at least to that 
precision. Furthermore, if one intends to achieve adia- 
baticity by increasing the switch-on time of the gate, so 
that / Jdt = (2n + X)-kTi with n > 1 pf, the error in 
J due to gate inaccuracies would accumulate in the in- 
tegral. It is thus crucial to precisely control the surface 
gates in order to minimize the possible errors. 

Placing the donors and therefore their electrons in the 
middle of a symmetric Si quantum well will help reduce 
the fast phase oscillations because only two of the six bulk 
Si conduction band valleys would contribute to the donor 
electron ground state [^), suppressing the magnitude of 
the interference effect. In a symmetric quantum well the 
splitting between the ground state and the next excited 
state (in analogy to the A\ and T donor electron states 
in bulk) is much smaller than in the bulk (crudely esti- 
mated to be about 3 meV from symmetry arguments, in 
contrast to nearly 12 meV in bulk Si; the smaller splitting 
implies that the adiabatic condition is somewhat harder 
to satisfy in the quantum well). Strains, electric fields, 
and/or asymmetry in the quantum well might help elim- 
inate the remaining valley degeneracy || so that all the 
fast oscillating factors in the exchange constant J may be 



removed, thus leading to a slowly varying J and a much 
easier control of the exchange gate. However, strains 
themselves are hardly controllable, while asymmetry in 
a quantum well may introduce additional complications 
that lead to decoherencc. Further investigations of these 
factors are still ongoing and the results will be reported 
elsewhere p"7[ . 

As we have demonstrated above, moving one of the 
donors to its second nearest neighbor position causes 
strong suppression in exchange coupling between the two 
donor electrons. This cancelation of exchange due to val- 
ley interference might be useful for isolating the neigh- 
boring qubits. In a Si QC, for most of the time the 
qubits should not "talk" to each other. Unwanted inter- 
action would lead to leak of information and decoherence. 
Therefore, the positions with vanishing exchange inter- 
action might provide a quieter environment compared to 
an arbitrary pair of positions for qubits. 
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